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Charged impurities in semiconductor quantum dots comprise one of the main obstacles to achiev- 
ing scalable fabrication and manipulation of singlet-triplet spin qubits. We theoretically show that 
using dots that contain several electrons each can help to overcome this problem through the screen- 
ing of the rough and noisy impurity potential by the excess electrons. We demonstrate how the 
desired screening properties turn on as the number of electrons is increased, and we characterize 
the properties of a double quantum dot singlet-triplet qubit for small odd numbers of electrons per 
dot. We show that the sensitivity of the multi-electron qubit to charge noise may be an order of 
magnitude smaller than that of the two-electron qubit. 



One of the most promising paths to scalable quantum 
computation is to use laterally defined double quantum 
dots (DQDs) in semiconductor heterostructures. The 
qubit is formed by the spin states of the two-electron 
DQD with total spin projection zero along the z axisJ^ 
Such a qubit is insensitive to spatially uniform magnetic 
field fluctuations, and, most importantly, amenable to 
fast electrical manipulation^^i^ Recent experiments have 
made tremendous advances along these lines, demon- 
strating single-qubit initialization, arbitrary manipula- 
tion, and single-shot readout, all within a fraction of the 
coherence time of the qubit. Preliminary steps toward 
an entangling two-qubit gate have also been reported.^ 
In principle, successful completion of that program leaves 
the (admittedly enormous) challenge of scaling up to 
large numbers of qubits as the last remaining hurdle in 
the fabrication of a practical quantum computer. 

However, a practical issue has emerged which threat- 
ens to be a crippling impediment to continued rapid 
progress. The semiconductor samples used to create 
quantum dots invariably contain a number of charge 
impurity centers, perhaps lO^'^cm"^ in GaAs systemsj^, 
Even if the charge on these centers can be frozen to avoid 
switching noise, their presence inhibits access to the one- 
electron-per-dot regime since the lowest energy states of 
the dot may be fragmented due to the roughened poten- 
tial landscape ii2r— This makes it difficult to find samples 
suitable for spin qubit realization. Furthermore, typically 
the impurities do introduce some switching noise^^^— 
so that even in samples in which the impurities are all 
far enough from the DQD that a two-electron singlet- 
triplet qubit can be accessed, the interdot exchange en- 
ergy is still subject to random fluctuation, leading to 
gate errors and decoherence.^®"'^° This necessitates op- 
erating in a parameter regime such that the sensitivity 
of the exchange energy to the charge noise is minimized, 
a so-called "sweet spot".^^ In general, the charge noise 
problem is even more pernicious when performing two- 
qubit operations directly mediated by the Coulomb in- 
teraction, and one must again seek a sweet spot.— 
However, in practice, this strategy may not be sufficient 
since one typically cannot optimize over all noise chan- 



nels simultaneously!^ 

An alternative approach to overcoming the charge 
noise problem is to define qubits with several electrons 
per dot, such that the "core" electrons (electrons paired 
up into singlets and filling up the lowest single-particle 
states comprising the multi-electron ground state) serve 
to screen out the charge impurities, while each dot con- 
tributes the spin of a lone "valence" electron to form 
the qubit. Experimental studies of spin blockade^i^ 
and excitation spectr have already been performed 
in multi-electron DQDs. On the theoretical side, early 
work analyzed some of the subtleties involved with qubit 
manipulation in the context of multi-electron dots,— 
highlighting, for example, the need for an external mag- 
netic field in the case of singlet-triplet qubits in order 
to ensure that the qubit is well defined. We shall see 
that the use of multi-electron dots for spin-based quan- 
tum computation merits further consideration because of 
the favorable behavior of such qubits in the presence of 
charge impurities. 

In this work, we demonstrate the effectiveness of mul- 
tiple electrons per dot in screening a single charged im- 
purity. We do this by performing a detailed numerical 
analysis employing the configuration interaction method 
with up to ten electrons in the single dot case and up to 
14 electrons in the case of DQDs. In the case of a single 
quantum dot, we compute the ground state energy twice: 
once for a multi-electron state which takes into account 
interactions between the electrons and the impurity, and 
a second time for a multi-electron state which ignores 
these interactions. This yields two ground state ener- 
gies whose difference constitutes a clear and convenient 
measure of the screening mechanism since it is directly 
related to the re-adjustment of the multi-electron wave- 
function in response to the impurity. We find that the 
effect of the screening grows quickly as the number of 
electrons in the dot is increased. 

In the case of DQDs, we calculate the exchange energy 
of a singlet-triplet qubit with three to seven electrons per 
dot and analyze its sensitivity to the charged impurity. 
To properly compare exchange energies obtained for dif- 
ferent numbers of electrons, we first define a valence elec- 
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tron tunneling rate, and as we change the number of elec- 
trons, we tune the DQD potential to keep this tunneling 
rate invariant. This helps to isolate true multi-electron 
physics from other effects which occur as a byproduct 
of adding electrons to the system. We find that as the 
number of electrons is increased, the exchange energy be- 
comes significantly less sensitive to the impurity. With 
five electrons per dot, for example, the sensitivity can be 
reduced by nearly an order of magnitude relative to the 
one-electron-per-dot case. Having five electrons per dot 
also appears substantially better than having three; this 
is due at least in part to the fact that, in a magnetic field 
(though not the high fields considered in Ref. [Sy) , there 
occurs some anomalous behavior arising from the spatial 
dependence of the phase of the relevant valence orbital.-^ 




FIG. 1: Model double well potential along the interdot axis. 
There is also a harmonic potential along y with frequency ujq. 
We have schematically shown 10 electrons filling the lowest 
Fock-Darwin orbitals in a finite magnetic field. 



In contrast to previous theoretical studies of multi- 
electron DQDs^"— or multi-electron qubitS j^^i^^'^^ our 
focus is on how charge noise affects singlet-triplet qubits 
in a DQD. However, our results also give some new in- 
sights pertaining to the multi-electron quantum dot sys- 
tem in the absence of charge impurities, and so has some 
direct bearing on these earlier works. In particular, some 
of these work o^^i^^ made use of the so-called "frozen- 
core" approximation, which is an approximation often 
employed in the context of numerical solutions (partic- 
ularly the configuration interaction method) to multi- 
electron quantum dot problems. In this approximation, 
one keeps the core electrons frozen in the lowest non- 
interacting states of the dot(s) and only allows config- 
urations where the valence electrons are free to occupy 
higher energy levels. In the course of our analysis of the 
charged impurity screening effect, we will implement a 
more sophisticated approximation scheme wherein exci- 
tations of the core electrons are also taken into account, 
allowing us to test the accuracy of the frozen-core ap- 
proximation. We find that the frozen-core approximation 
becomes significantly worse as the number of electrons is 
increased. Although the method we employ constitutes a 
marked improvement over previous approaches, we focus 
on qualitative results and general trends since the micro- 
scopic details of the actual experimental potential and 
noise are not precisely known. 



The structure of the paper is as follows: In Sec. HI we 
explain our numerical methods and examine the valid- 
ity of the approximations employed. In Sec. |TT1 we show 
the onset of screening effects in a multi-electron single 
quantum dot. In Sec. IIIIl we characterize multi-electron 
singlet-triplet qubits in DQDs and their sensitivity to 
charge noise. Finally, we conclude in Sec. IIVI In Ap- 
pendix |X1 we review Fock-Darwin states to keep the pa- 
per self-contained and to fix conventions. Appendices iBl 
and [Cl contain technical details that facilitate the numer- 
ical computations of the multi-electron spectra. 



I. NUMERICAL METHODS 

A. Basis states 

All of our results are obtained in a configuration- 
interaction (CI) approach using molecular orbitals. The 
Hamiltonian for N electrons in confining potential V (r) , 
in the presence of M static charge impurities located at 
positions R^-, and with a perpendicular magnetic field 
B = Bz = V X A, is 
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where m* is the effective mass of the electrons, Zi^p is 
the impurity strength, e is the dielectric constant of the 
semiconductor, g* is the effective electron g-factor, /is 
is the Bohr magneton, and S is the total electronic spin. 
In this work we will take parameters relevant to GaAs 
systems: m*=0.067me, g* = -0.44, e = 13.1. We neglect 
perturbative corrections to the Hamiltonian such as spin- 
orbit coupling. When including the effects of an impurity, 
we will consistently take Zi^p =^ 1 . 

We will assume tight confinement along the z-direction 
in the depletion layer of a heterostructure interface and 
consider the resultant two-dimensional electron gas. For 
simplicity, we will use a (bi-) quadratic in-plane poten- 
tial to model the lateral gate-defined confinement of a 
(double) quantum dot. 



V {x, y) — —m*ujQMin{{x ■ 



(2) 

as shown in Fig. [1] We have checked that, as expected, 
the discontinuity in the derivative of the potential does 
not have any qualitative effect. Although the actual con- 
finement potential depends on the experimental details 
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of the gate geometry and voltage, the biquadratic ap- 
proximation has been justified by comparing to an ex- 
act numerical solution of Poisson's equation for a realis- 
tic setup. The appropriate single-particle orbitals are 
then the Fock-Darwin states centered at the minima of 
each well, r = ^xoi., labeled with a ± sign, respectively, 
and carrying principal quantum number n = 0, 1, 2, ... and 
magnetic quantum number ra = —n, —n + 2, n — 2,n, 
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where £„ = £b/ (1/4 + ^2/^2)1/4^ ^ ^l^/Zs, = 
eB/ 171*0, (x) is an associated Laguerre polyno- 
mial, and we have taken the symmetric gauge, A = 
B/2 (— yx -I- xy). The associated single-particle energies 
are 



En,m = (n+l) 




(4) 



The derivation of the Fock-Darwin states and spectrum is 
reviewed in Appendix [X] In analogy with the hydrogen 
atom, we refer to the {n,m) = (0,0) state as the "5"' 
orbital, the {n,m) = (1,±1) state as the "P±" orbital, 
and so on. We will similarly refer to groups of states 
with the same n as "shells." 

The basis states used to expand the many-electron 
wavefunction are Slater determinants formed from the 
single- particle orbitals, 

(5) 

where a is the electronic spin variable, x (o") is the associ- 
ated spin wavefunction, and V is the antisymmetrization 
operator. See the appendices for a complete presentation 
of the basis states and useful matrix elements between 
them. 



B. Truncation and convergence 

While the exact ground state can be built from the 
complete set of configurations using all combinations of 
Fock-Darwin single-particle states, obviously in practice 
one must truncate the space of possible configurations. 
The analogy with the hydrogen atom spectrum naturally 
leads one to the idea that perhaps, in dealing with multi- 
particle states, one can treat the core electrons as be- 
ing effectively inert. In other words, one can consider 
the "frozen-core" approximation in which one keeps only 
multi-particle states for which all but two of the electrons 



fill up the Fock-Darwin states below the valence orbital.— 
Another technique is to simply neglect all configurations 
involving single-particle orbitals above some cutoff level. 
This is referred to as "full CI" within the orbital cut- 
off. A more efficient method is to neglect all configu- 
rations whose non-interacting energy is higher than the 
non-interacting ground state by some cutoff energy. We 
have used each of these approximations in our numeri- 
cal calculations for comparison. Particularly for several 
electrons and in the presence of an impurity, the frozen- 
core approximation may not allow enough flexibility in 
the multi-electron wavefunction to accurately represent 
the ground state. 

In a single dot with typical radius on the order of 
100 irni^"— the on-site Coulomb energy is much greater 
than the non-interacting level spacing, so using the basis 
built from non-interacting Fock-Darwin orbitals requires 
a very large number of excited configurations to be kept 
in the CI calculation. Since the ratio of Coulomb energy 
to level spacing is proportional to dot size, we take a rel- 
atively tight confinement potential with oscillator length 
between 6 — 10 nni to aid convergence. Even for these 
strongly confined electrons, though, the interaction en- 
ergy scale is comparable to the confinement energy scale 
and we must keep many excited configurations. Exper- 
imental dot radii are usually considerably larger^"— but 
we expect that the features will be qualitatively similar. 

In Fig. [21 we show the convergence of the ground 
state energy for three and five electrons in the absence 
of impurities within a full CI calculation as the number 
of shells kept increases. For five electrons, keeping six 
shells amounts to keeping nearly 20,000 configurations. 
We have chosen values of the magnetic field such that 
LOc — 2uj() I -\/99, which, although it corresponds to a ~ 3T 
field for our unusually small dots, for a typical dot size 
corresponds to a field of only tens or hundreds of mT. 
(The convergence results are specific to the choice of pa- 
rameters. The particular set of parameters taken here is 
in analogy to the DQD potential discussed in Sec. lIIIBl ) 
Also shown are the energies obtained using a frozen-core 
approximation in which only the outermost electron is 
free to occupy the excited shells. Clearly the frozen-core 
results are significantly improved upon by allowing con- 
figurations with excited core electrons. Again, this is due 
to the strong electron-electron interactions. 

In Fig. 121 we keep all configurations with non- 
interacting energy within a given cutoff excitation energy, 
-Ec, from the non-interacting ground state energy. For ex- 
ample, when i5c = 0, only the ground-state configuration 
is kept. We show the convergence of the ground state 
energy versus cutoff for three and five electrons. We use 
hijJc/2 as a natural unit because that is the smallest incre- 
ment by which the single-particle energies of Eq. (U) can 
change. The step-like behavior is due to the discreteness 
of the spectrum and the fact that angular momentum is a 
good quantum number: Increasing the cutoff energy only 
increases the basis size at a discrete set of points, and a 
new configuration with the correct angular momentum 
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FIG. 2: (Color online.) Ground state energy vs. number of 
shells kept in the full CI calculation for three (a) and five (b) 
electrons in a harmonic trap with frequency fejo^ 24.26 meV 
(a) and /laio = 21.79 meV (b) and perpendicular magnetic field 
B = 2.8 T (a) and B = 2.5T (b) such that a;c = 2a;o/y99. 
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FIG. 3: (Color online.) Ground state energy vs. cutoff in the 
energetically truncated CI calculation for three (a) and five 
(b) electrons in a harmonic trap with parameters as in Fig. [2] 



further in Sec. UlTC] 



is added to the basis only when the cutoff permits con- 
figurations with an even increment in the total principal 
quantum number, n,;. For the parameters chosen, this 
occurs every 20hu!c/2. We see that taking an excitation 
cutoff of Ec ~ 80fta;c/2 gives the same energy as a full 
Cf approach with six shells, while only keeping about a 
tenth of the number of configurations. However, even a 
cutoff of Ec ^ 20hajc/2 already captures the dominant 
corrections to the non-interacting state. Somewhat sim- 
ilar convergence properties were noted in Ref. [s^ 

Similar calculations for the tunnel-coupled DQD sys- 
tem (with twice as many electrons) become computation- 
ally demanding at high cutoffs, but from the above results 
for a single dot (or equivalently, for well separated DQDs) 
we expect reasonable convergence is attained already for 
full CI up to the F orbitals or for a cutoff of Ec ~ 2Qhajc/2. 
Since full CI up to F orbitals quickly becomes unwieldy 
for several electrons, we shall primarily use the cutoff ap- 
proach for the DQD calculations. For five electrons per 
dot, taking numerically tractable cutoffs on the order of 
20huJc/2 corresponds to keeping roughly the lowest thou- 
sand non-interacting DQD configurations, which we note 
to be a significant improvement on the common approx- 
imation of keeping only S", P, and sometimes D shells 
with a frozen core. We will discuss DQD convergence 



II. ONSET OF SCREENING FOR 
MULTI-ELECTRON DOTS 

We now consider the multi-electron single quantum 
dot in the presence of a single charged impurity. Upon 
turning on the impurity potential, the electrons redis- 
tribute themselves to minimize their total energy. The 
amount by which the energy in the presence of the impu- 
rity changes due the rearrangement of the multi-electron 
wavefunction is 

A= (*|i/|*)-('J'o|ff|'J'o), (6) 

where \^) and |^o) are the ground states with and 
without an impurity present, respectively, and H is the 
Hamiltonian with an impurity. This quantity is plotted 
in Figs, m and [S] for different positions of the impurity. 
The calculation is performed with cutoffs in the range 
Ec — A^1tWc/2 to Ec = llOhuJc/2, with lower cutoffs used 
for higher electron numbers. For Fig. 21 the minima of 
the combined impurity plus harmonic potential remains 
at the origin, and A may be thought of as a measure of 
the change in the shape of the wavefunction. For Fig. [SJ 
the minima of the potential shifts, and A is much larger 
due to a trivial contribution coming from the translation 
of the wavefunction. 
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FIG. 4: Rearrangement energy A (see text) vs. number of 
electrons for an impurity on the z-axis. Here the impurity is 
repulsive, and is located 17 nm (top) or 34 nm (bottom) from 
the dot center. There is a perpendicular field B = 2.8T and 
the trap frequency is huo — 24:.26meV . 



FIG. 5: Rearrangement energy A (see text) vs. number of 
electrons for an impurity on the x-axis. Parameters are the 
same as in Fig. [l] 



III. MULTI-ELECTRON SINGLET-TRIPLET 
QUBIT IN THE PRESENCE OF CHARGE NOISE 



In both cases, the linear part of the curves may be un- 
derstood as arising from the impurity pushing the elec- 
trons away from the center of the harmonic potential. 
However, note that one begins to see marked deviations 
from linearity already at six or seven electronSi^. This 
implies that nontrivial screening effects are already set- 
ting in whereby the first few electrons form a screening 
core which rearranges to partially shield any additional 
outer electrons from the impurity potential. In other 
words, defining a multi-electron spin qubit using even a 
small number of electrons per dot is advantageous since 
screening reduces the effect of the impurity. This is con- 
sistent with previous qualitative predictions of screening 
effects in the optical absorption of few-electron dots'^* 
and in the high magnetic field limit of a six-electron dot 
with a repulsive impurity.— We shall see in Sec. lIII Dj that 
possible advantages are not limited to qubit fabrication 
in the presence of static impurities, but also apply to 
qubit operation in the presence of fluctuating impurities. 



The screening discussed in the previous section may aid 
the fabrication of multi-electron quantum dot systems in 
disordered background potential landscapes that would 
fragment or destroy single-electron dots. For the pur- 
poses of quantum computation, we now wish to discuss 
the properties of a multi-electron singlet-triplet qubit 
that can take advantage of this screening. 



A. Qubit definition 

Consider a DQD with an odd number of electrons 
in each dot. In the absence of impurities and without 
the Coulomb interaction between electrons, the electrons 
forming the ground state will pair up into singlets and 
inhabit the lowest Fock-Darwin states. If the tower of 
Fock-Darwin states in each dot is non-degenerate (which 
requires a non-zero magnetic field), then for odd num- 
bers of electrons in each dot it will always be the case 
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that there is exactly one highest occupied valence state 
in each dot, and the valence state in each dot contains 
only one "valence" electron. For = 2, 6, 10, 14 electrons 
(iV is the total number of electrons in both dots), the 
valence state is the S {ti = 0, m — 0), P_ {n = l, m — ~l), 
P+ (ri = l, m = l), and D- {n = 2, m = — 2) orbital, re- 
spectively. The spins of the valence electrons can then be 
used to define a singlet-triplet qubit in exact analogy to 
the standard two-electron casci^ Even when the intuitive 
picture of an inert core breaks down, so that there is no 
sharp distinction between valence and core electrons, the 
qubit is still defined in terms of the lowest two many-body 
energy levels and the low-energy spectrum of the DQD 
retains a good qubit subspace as will be shown below. 



B. Parameter choices 

As mentioned previously, for better convergence we 
take small dots with characteristic radius between 6—10 
nm. We choose the inter-dot distance to be 40 nm in 
order to have distinct dots that still have appreciable 
tunnel coupling and exchange. This is a smaller inter- 
dot distance than is typical for experiments^^— but that 
is simply because the dots themselves are smaller here. 
The qualitative behavior should be similar. 

When comparing cases with different electron num- 
bers, one should first isolate the intrinsic multi-electron 
effects from trivial filling effects. For example, for a 
given double well potential, a six-electron singlet-triplet 
qubit will trivially have a larger exchange energy than a 
two-electron one simply because the valence electrons lie 
higher in the well and thus see a lower barrier to tunnel- 
ing. This hardly makes for an enlightening comparison 
between the two- and six-electron cases, since the domi- 
nant effect of naively adding more electrons is effectively 
nothing more than a change in the central barrier height. 
We will instead compare different fillings on a more nu- 
anced and consistent basis by adjusting the magnetic field 
and the confinement potential such that the tunneling be- 
tween valence shells and the ratio of orbital level spacings 
remain constant as the filling is varied. 

More precisely, we define the valence tunnel coupling 
for N electrons in two dots to be the energy splitting 
in the frozen-core approximation between the two lowest 
states for A^— 1 interacting electrons in two dots, the idea 
being that this energy splitting characterizes the tun- 
neling of an electron in the valence state of one dot to 
the valence state of the other. This tunneling incorpo- 
rates the Coulomb interactions with all core electrons, 
but does not include interactions with impurities. By ad- 
justing the confinement strength and the magnetic field, 
we fixed the value of this tunneling to be 0.2 meV for 
all cases. This particular value was chosen to achieve 
good numerical convergence while remaining in the win- 
dow of experimentally relevant exchange energies, which 
are typically on the order of a few fieYi^ in our setup, this 
corresponds to tunneling rates on the order of hundreds 
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TABLE I: Energies of low-lying Fock-Darwin shells for 
ojc/ojo = 2/^99. 



of fieV. 

Before quoting the parameters needed to obtain this 
tunneling for the cases of7V = 2,6,10,14 electrons, we 
first need to say a few words about the spectrum. We 
mentioned above that we need an external magnetic field 
to obtain a non-degenerate spectrum in each dot. How- 
ever, we did not specify precisely what value this i?-field 
should take. We want the i?-field to be large enough to 
avoid possible near-degeneracy effects, but at the same 
time, we do not want it too large because then different 
shells would start to overlap in energy. Of course, the 
latter effect happens for any finite i3-field at sufficiently 
large n. In practice though, very high shells should not 
be relevant, and for the most part we will only keep states 
up through the n = 3 shell. (For some of the = 14 re- 
sults, we will find it necessary to also include some of the 
n = 4 states, but this has little bearing on the present 
considerations and will be further clarified below in the 
context of the particular results in question.) It then suf- 
fices to make the B-field small enough so as to avoid any 
overlap between the n = 3 "F" shell and the n = 4 "G" 
shell. We can satisfy all the above criteria by choosing 
the cyclotron frequency 



leading to the Fock-Darwin spectrum of Table U and non- 
interacting mult i- particle energies 

N N . 

^{Tiita = E = E (10"' + 10 + m,) 

1=1 j=i 

(8) 

For this choice of parameters, the spacings between states 
of the same shell are smaller than the spacings between 
shells, but not drastically so. Also notice that we are 
tying the _B-field to the size of the wells. This means 
that when we change the number of electrons, we will 
simultaneously change both the B-field and well size in 
such a way that the valence tunneling (as defined above) 
remains constant at 0.2 meV. 

With these considerations, we calculated the frozen- 
core spectrum of 1, 5, 9, and 13 electrons in a double 
well with interdot spacing of 40 nm to obtain the charac- 
teristic tunneling energies shown in Fig. [5] as a function of 
confinement strength. Note that the unphysical behavior 
of the tunneling for small ujq is an artifact of the trunca- 
tion of the configuration space, an approximation which 
clearly breaks down for weak confinement. From Fig. [S] 
we obtain the parameters (listed in Table |n| necessary 
to equalize the valence tunneling at 0.2 meV for the 2, 
6, 10, and 14-electron singlet-triplet qubits we wish to 
consider. 
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FIG. 6: (Color online.) Tunnel coupling versus well depth 
for 1, 5, 9, and 13 electrons in two dots in the frozen-core 
approximation. 
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TABLE II: Parameters for which the valence tunneling is t = 
0.2 meV with interdot distance fixed at 40 nm. 

Note also from Fig. |6] that the valence tunneling is 
actually smaller for nine electrons, where the valence or- 
bital is P+, than for five electrons, where the valence 
orbital is P_. It is rather unusual and counterintu- 
itive that a higher energy orbital has a smaller tunneling 
rate! We interpret this as a consequence of the phase 
^±ixoy/i2eg) appearing in the n = 1 shell hopping- type 
integrals, / dr (^dx(j)t ±i) dx4'i±i, "^^^ the magnetic 
field (see Eq. ([3|)). This allows contributions from the 
prefactor in the integrand which is an odd function of y. 
Since, for xq > £o, the odd term has a -|- (— ) sign for 
m<0 (m>0), the magnetic field diminishes the tunnel- 
ing of the higher P orbital and enhances that of the lower 
one. Somewhat similar considerations have been noted 
in recent experimental worki^ 



C. Qubit characterization 

We characterize these multi-electron singlet-triplet 
qubits by calculating the energy splitting, J, between 
the two lowest levels in the absence of a charged impu- 
rity. (These are well separated from the other levels, as 
shown below, and form a good qubit subspace.) Previ- 
ous studies considered this splitting, which is an effec- 
tive exchange energy, for six electrons using a frozen- 
core approximation^SS, We will improve upon these early 
results both by increasing the number of electrons to in- 
clude the = 10 and TV = 14 cases, and by relaxing 
the assumption that the core electrons are frozen. In 
order to both test and improve upon the frozen-core ap- 
proximation, we will use the cutoff approximation as in 



Sec. lIBl keeping all multi-particle states that have a non- 
interacting energy below a certain cutoff excitation en- 
ergy, Ec, from the non-interacting ground state. When 
Ec = 0, we obviously have only the four multi-particle 
states which correspond to the four ways of distributing 
the two valence electrons among the two valence orbitals. 
As the cutoff is increased, the number of states included 
quickly increases at a rate which grows with N. 

In Fig. [7] we show the convergence of the lowest two 
eigenenergies, and in Fig. [5] we show their difference, J. 
In the two-electron case it appears that J is well con- 
verged for cutoffs larger than Ec — 20huJc/2, in agree- 
ment with our expectations from the single-dot results 
of Sec. II Bl It is also apparent in Figs. [7] and H] that a 
substantial jump occurs at Ec = 2QhiOc/'^ for all num- 
bers of electrons. The significance of this special cut- 
off value (which depends on our choice of parameters) 
is that it corresponds to the point at which the lowest 
excited state with the same magnetic quantum number 
TO as the valence orbital is included in the cutoff scheme. 
When the valence orbital has quantum numbers (n,TO), 
this newly added state has quantum numbers {n + 2, to). 
Coulomb matrix elements involving this state and the 
valence orbital are quite large compared with matrix el- 
ements between states with different values of m. The 
reason can be seen by considering a single dot (or equiv- 
alently, a DQD with infinite interwell separation). In 
that case, angular momentum is a good quantum num- 
ber. Thus the ground state must be built from a set 
of orbitals which all have the same angular momentum. 
Thus, if one considers the convergence of ground state 
energy as orbitals are added to the basis set one-by-one, 
the energy should only be lowered when an orbital with 
magnetic quantum number to is added to the basis. Now, 
angular momentum is no longer conserved for a coupled 
DQD, and the ground state contains contributions from 
all orbitals. However, for the physically relevant interdot 
distance we have taken, the ground state still favors the 
set of orbitals discussed above. Hence the energy changes 
more when we add the (n + 2, to) orbital than when we 
add others. 

These observations highlight the importance of using 
cutoff energies above Ec = 20huJc/2. As Ec is increased, it 
eventually becomes large enough that Fock-Darwin states 
beyond those kept in the CI calculations become excited, 
and it would be inconsistent to increase Ec any further 
since we would be adding higher-energy configurations 
without first including all configurations with lower en- 
ergy. For the cases = 2, 6, 10, it suffices to keep up 
through the F shell in order to reach the Ec > 20hujc/2 
regime, and the upper bound on Ec arising from these 
considerations is fixed by the energy difference between 

the valence state and the G orbital (and so depends on 

the number of electrons). For the iV=:14 case, however, 
it is necessary to include two orbitals from the G shell in 
order to achieve Ec > 20huJc/2. Thus for iV = 2, 6, 10, 14 
we must have Ec/ {huJc/2) < 36,27,25,22, respectively. 
Although it is evident from Figs. [7] and [5] that small 
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fluctuations in the lowest energies and J continue to 
arise beyond = 20fta;c/2, there are indications that 
no additional large jumps occur as Ec is increased fur- 
ther. In particular, we have noticed that such jumps 
are also apparent in the single-particle matrix elements. 
The angular-momentum considerations above would sug- 
gest that such jumps would occur at integer multiples 
of 20hijJc/2; however, we have computed single-particle 
matrix elements including up through the n = 6 shell 
and found that no further jumps arise. Although a sim- 
ilar calculation for multi-particle matrix elements would 
be computationally very challenging, these single-particle 
results do at least suggest that the curves in Figs. [7] and 
m are reasonably well converged. 

We have used the parameters from the previous sub- 
section so we can directly compare the exchange ener- 
gies for different electron numbers at a fixed tunneling 
rate, t ~ 0.2 meV. We see from Fig. [5] that, in con- 
trast to the zero magnetic field case within a frozen-core 
approximation;^ the exchange energy generally increases 
for larger numbers of electrons. 

Fig. ini shows how the low-lying spectrum behaves as a 
function of harmonic well frequency, wq, with magnetic 
field chosen such that uJc/ujq = 2/y/99 is held constant. 
We have subtracted off a trivial linear dependence on ujq 
arising from the non-interacting ground state energy us- 
ing Eq. ([5]). Only the lowest states are shown; in reality 
there are of course many closely-spaced excited states in 
the empty upper portion of the plots. Note that, as ex- 
pected, the lowest pair of states (indistinguishable on the 
scale shown) are isolated from the other states by a large 
energy splitting due to the tight confinement and strong 
interactions. The singlet-triplet splitting, J, is shown 
in the insets, including results using the frozen-core ap- 
proximation for comparison. Generally, the frozen-core 
approximation is valid in the limit of very large trap fre- 
quencies such that the Coulomb interactions become neg- 
ligible. However, for realistic dot sizes, the frozen-core re- 
striction underestimates the exchange considerably. The 
discrepancy becomes more pronounced as the dot filling 
is increased. 



D. Qubit operation 

Finally, in addition to the benefits in fabrication, we 
now show that the screening is likewise beneficial in 
the operation of singlet-triplet qubits in the presence of 
charge noise. Random telegraph fluctuations in the occu- 
pation of a charge impurity center results in fluctuations, 
6 J, in the exchange energy, inducing singlet-triplet deco- 
herence and gate errors. Recent experiments have made 
remarkable progress in countering dephasing due to hy- 
perfine coupling to the nuclear spin bath, extending the 
decoherence time to ~ lOO^s,— leaving charge fluctua- 
tions as the dominant source of error. This presents a 
formidable obstacle to current efforts to perform high- 
fidelity logical gates. Although one may be able to find 
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FIG. 7: (Color online.) Ground and first excited state ener- 
gies (left) for 2, 6, 10, and 14 (top to bottom) electrons and 
the difference between the two (right) vs. energy cutoff. 



"sweet spots" where 5 J is minimized for a given noise 
channel^ e.g., fluctuation in interwell detuning, one 
generally cannot simultaneously protect against fluctu- 
ations in the other model parameters.'^- Furthermore, in 
a large-scale system, the fine-tuning required to ensure 
that the whole ensemble is operating at a sweet spot 
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FIG. 8: (Color online.) Splitting between ground and first 
excited state energies vs. energy cutolf for 2, 6, 10, and 14 
electrons. 



would seem to be prohibitive. On the other hand, using 
muhi-electron DQDs in order to reduce 5 J via screening 
is a relatively general and straightforward approach. It is 
thus desirable to consider the effectiveness of the screen- 
ing in protecting against noise from charge fluctuators. 

Decoherence from the charge noise is essentially deter- 
mined by the average switching time, r, for all but the 
most minute values of 5 J on the order of h/r. Gate er- 
rors, though, are dominated by the relative fluctuation in 
the exchange energy, 5J/J. We calculate this quantity 
for a single repulsive charge impurity center located di- 
rectly above one of the dots, neglecting any image charge. 
This is only meant to give a qualitative picture of gate 
errors due to charge noise. Although it would be no more 
difficult to perform the calculation with any given impu- 
rity potential, in the absence of direct knowledge of the 
impurity positions and switching times relevant to a spe- 
cific sample, any choice is purely arbitrary. As shown in 
Fig. [TUl the relative fluctuation in J for a multi-electron 
qubit may be reduced by as much as an order of magni- 
tude compared to the single-electron-per-dot case, which 
would result in a dramatic reduction of the charge noise 
induced infidelity. Figure [TU] shows the results of using 
the best energy cutoff approximation in each case. 

Generally speaking, more electrons correspond to less 
sensitivity to the charge fluctuator. In the ten-electron 
case, though, the sensitivity appears anomalously small. 
This may be due to the valence electrons residing in the 
P+ orbitals, which have a suppressed tunneling rate, as 
discussed above. As a result of this orbital effect, the 
trapping frequency, wq, does not have to be increased as 
much to match the two-electron S orbital tunneling rate. 
This in turn implies that the core electrons have more 
freedom to screen the impurity. 

Fig. [TT] shows the convergence of 5 J as the energy cut- 
off Ec approaches its maximal value, as well as full CI 
results for the two-electron case and frozen-core results 
for the other cases. It is not guaranteed that our results 
for the exchange energy are fully converged, especially 
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FIG. 10: (Color online.) Fractional change in exchange energy 
due to impurity vs. impurity distance for different numbers 
of electrons using the energy cutoff approximation with — 
35, 26, 24, 17 for iV = 2, 6, 10, 14 electrons respectively. 



for larger numbers of electrons, but Fig. II II demonstrates 
that the qualitative trends apparent in Fig. [10] are re- 
liable, most notably the decrease in sensitivity to the 
charged impurity for larger numbers of electrons. These 
figures also illustrate the worsening of the frozen-core ap- 
proximation as electrons are added to the system. 



IV. CONCLUSIONS 

The presence of randomly positioned and temporally 
fluctuating charge impurity centers in even the cleanest 
semiconductor samples is currently the most prominent 
roadblock to experimental realization of precisely con- 
trolled, long-lived semiconductor spin qubits. It is an 
open question to what extent and in what direction the 
path to scalable spin quantum computation will be al- 
tered by this roadblock. On one hand, one could pur- 
sue ultraclean samples with impurity concentration un- 
der lO^^cm"'^.^ On the other hand, one could seek to 
engineer qubits that are less sensitive to the presence of 
charge impurities, as we have discussed in this paper. 
However, since the standard two-qubit gate for singlet- 
triplet qubits relies on the Coulomb interaction;^- engi- 
neering a qubit that is less sensitive to fluctuations in 
the electrostatic environment also diminishes the abil- 
ity to manipulate the many-qubit state via the standard 
techniques. This is, of course, a familiar conundrum in 
all proposals of quantum computation. The outlook is 
nevertheless quite promising, since random fluctuations 
in the electrostatic potential can be suppressed by bias 
cooling^'^ or by adding a negatively biased insulated top 
gatei^ in order to suppress leakage of electrons from the 
gates into the two-dimensional electron gas. Further- 
more, intentional changes in the electrostatic potential 
for the purposes of two-qubit operations can be enhanced 
by a floating interdot capacitive gate^Mi It is therefore 
reasonable to sacrifice a little of the abundance of con- 



trollability in order to gain a measure of immunity to the 
impurity background. 

In this paper, we have explored the ramifications of 
using multi-electron quantum dots in order to diminish 
the effect of the random impurity potential. The numeri- 
cal calculations in this work provide qualitative guidance 
to experiments regarding the characteristics of a DQD 
loaded with 2, 6, 10, or 14 electrons and how its behavior 
depends on dot size and nearby charge impurity centers. 
Our method of approximation consistently accounts for 
deviations of the many-body wavefunction from an effec- 
tive two-electron, frozen-core treatment, and our qualita- 
tive observations are shown to be robust and physically 
sensible. Pursuing quantitatively more precise results 
would be fairly meaningless since we are using a sim- 
plified model of the confinement potential. A precise cal- 
culation would require a fully self-consistent Schrodinger- 
Poisson solution for a particular gate geometry and a spe- 
cific impurity distribution42ii^ However, the simpler and 
more general model potential has previously been found 
to be a good approximation to the exact one^^ In any 
case, the essential screening physics discussed above is a 
general feature which should not depend sensitively on 
the exact form of the confinement. A more important is- 
sue is the fact that our calculations were performed in the 
case of an unusually tight confinement in order to reduce 
the relative strength of the Coulomb interactions and aid 
convergence. This precludes a direct connection to recent 
experiments using dots several times larger. However, for 
larger dots the features we have discussed should become 
even more pronounced as the deviations from the frozen- 
core approximation become more important. 

In this work, we have characterized multi-electron 
qubits in the presence of a charged impurity. We have 
not discussed details of initialization, manipulation, and 
readout. We envision these being performed as in the 
two-electron cas©5. by tilting the double-well so that the 
valence electrons begin to shift onto the same site, result- 
ing in larger overlap and exchange. These are separate 
issues which, although not expected to be problematic, 
might require some care and could constitute the sub- 
ject of future investigation. We have demonstrated that 
using multi-electron DQDs in a finite magnetic field to 
form singlet-triplet qubits is a viable option that does not 
suffer from the suppression of the exchange energy that 
one might expect based on previous results in a more 
restricted approximation^^S Not only are multi-electron 
qubits viable, they are also an attractive option due to 
nascent screening of the rough background by the few 
"core" electrons in each dot. This results in reduced 
sensitivity to random static impurity potentials, allow- 
ing easier fabrication, as well as reduced sensitivity to 
fluctuations in the impurity potential, facilitating more 
reliable single-qubit manipulation. 

This work is supported by LPS-NSA and lARPA. 
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FIG. 11: (Color online.) Fractional change in exchange energy due to impurity vs. impurity distance for (a) 2, (b) 6, (c) 10, 
and (d) 14 electrons with various energy cutoffs. 



Appendix A: Fock-Darwin states 
1. One harmonic well 

The Fock-Darwin Hamiltonian has the form 



H = 



^{-m+lA)\lm*.y, (Al) 



where r = ^ -\- 'ip' is the 2d radial coordinate, m* is 
the effective electron mass, and the vector potential is 
given by 



A = -— yx+ — xy. 

Introducing the cyclotron frequency, 

eB 
m*c 

we can write the Hamiltonian as 



(A2) 



(A3) 



H 



ft2 



1 



2^.-V^ + ^™*K+-cV4).^ + ^f4, (A4) 
where is the z-component of the angular momentum: 
4 = -ih{xdy -yd.,). (A5) 



It will be useful to consider the Fock-Darwin Hamiltonian 
in dimensionless complex coordinates, defined by: 



X + iy _ X — iy 
z = — — — , z — 



V2£o 



where the length scale io is 

and the magnetic length is given by 




(A7) 



(A8) 



Defining the following set of creation and annihilation 
operators, 

^{z + d), at = --i=(z-a), 

b=l=iz + d), b^^^iz-d), (A9) 
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which have the commutation relations 

[a,at] = l, [6,6t]^l, (AlO) 

with aU other commutators vanishing, the Hamiltonian 
becomes 

H = huj+{a^a + 1/2) + hu-ib'' b + 1/2), (All) 



with 



(A12) 



The ground state of this HamiUonian hves in the kernel 
of both a and b and so has the (normalized) wavefunction 



(A13) 



^0 = \/ -e 



It is trivial to check that atpo = bipQ = 0. The commuta- 
tion relations (|A10p imply that 

b'^b{b^T = n(6^)" + (fo^)"+^6, (A14) 
so that the class of (normalized) functions 



l==(at)"+(6t)"-V'o (A15) 



+ !n_! 



are eigenfunctions of the Hamiltonian with eigenvalues 

En^,n^ = haj+{n+ + 1/2) + fiw_(n_ + 1/2). (A16) 

The quantum numbers n+ and n_ are nonnegative inte- 
gers. If we define 



n = 71+ + n_ 



(A17) 



then n is a nonnegative integer, and m takes values from 
— n to n in steps of 2. The spectrum in terms of n (prin- 
cipal quantum number) and m (magnetic quantum num- 
ber) is 



E„ 



{n + l)h^uj^ +uj'^/4 + mnujc/2. (A18) 



this coordinate shift simply by defining new complex co- 
ordinates: 



z± = X ± Xq + iy, z± = X ± xq — iy. 



(A19) 



However, this redefinition moves not only the location of 
the quantum dot but also changes the vector potential 
A, Eq. (jA2p . We can correct this by also performing a 
gauge transformation on A: 



(A20) 



At 2"2^oy, 



which enables us to keep A fixed to the form in Eq. (jA2j) 
while we move the harmonic well in the x-direction. Of 
course, a gauge transformation affects not only the vector 
potential but also the wavefunction of the electron, and 
in particular does so in such a way that the wavefunction 
picks up a phase: 



(A21) 



The phase (p can be fixed by going back to Eq. (|Aip . im- 
plementing the gauge transformation on A, and demand- 
ing that the new wavefunction is still an eigenfunction of 
the Hamiltonian, leading to the condition 



This condition has the solution 



eBxQ 



(A22) 



(A23) 



where the ± corresponds to the dot at a; = T^o, or in 
other words, the right dot has the minus sign in the phase, 
the left dot has the plus. Notice that this phase is the 
same for all states living in the same dot. The bottom 
line is that we may continue to use the same Fock-Darwin 
states as before, but with z and z now defined according 
to Eq. (|A19p and with the overall phase factor e"^± in- 
cluded, where ip± is given in Eq. (|X23l) . Also note that, 
when using the raising operators to generate states, the 
phase is to be included after the operators have been ap- 
plied to the phaseless 4'o'- 



2. Two wells located at {x,y) = (±a;o,0) 

In the case of two quantum dots, we need to modify 
slightly the Fock-Darwin states found above. Obviously, 
the above results are valid no matter where the dot is 
located so long as our coordinates are defined with re- 
spect to the center of the dot. In the case of two lat- 
eral quantum dots, however, we want to define the co- 
ordinates with respect to the point halfway in between 
the dots. In these coordinates, the dots are located at 
{x,y) = {±xo,0). One might be tempted to implement 



V2 



^(al)"+(6L)"-e 



(A24) 



Here, the -l-(— ) index gives a state in the left (right) dot. 
We have taken the liberty of defining the operators a^. 
and bj., which are just the usual creation operators but 
with z and z replaced by z± and z±. In terms of n and 
m quantum numbers, we have 



Xx,y). 



(A25) 
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Appendix B: One and two-particle matrix elements 
1. Integrals 

Matrix elements involving products of Fock-Darwin 
states can be simplified by exploiting the fact that all 
such states arc generated from e^^^. (In this section, we 
suppress the ± indices indicating to which dot the state 
belongs.) Since and are linear differential operators 
acting on e~^^, the Fock-Darwin wavefunctions have the 



form P(z, z)e~^^ where P{z, z) is some bivariate polyno- 
mial. This in turn means that we can also generate the 
states by instead starting from the generator 

(l){z, z, c, d) = ^^6"^^+'=^+''^ (Bl) 

and differentiating with respect to c and d. The explicit 
form of the differential operator which yields the Fock- 
Darwin states is 



n"+'"- — 



/g:(_i)"-2("+-"-)/2ar-"- Efcio (1) 



(-2)'' 



(n+- 



I — 



(-2)'° 



' (n_-n++fe)!^d^c ' 



(B2) 



and in terms of this operator, the states are 



V'n+,n-(^,-2) = (l){z,Z,C,d) 



d=—a 



with 



ip± eBxo 

Q! = -Z— = ±— TT — • 

2y Ahc 



(B3) 



(B4) 



(On the right-hand side of this expression, c is speed of 
light and not the parameter appearing in (j){z,z).) Now 
consider an arbitrary single-particle matrix element com- 
prised of some operator A sandwiched between two Fock- 
Darwin states: 



{n+,n-\A\n'+,n'_) = 11 dxdyipl^^^_{z, z)Atpn'_^y_{z, z) = D"+'" // dxdy(p{z, z,c,d)A(j){z, z,c' ,d') 



c=a ' 
d=—a 

(}'--<y' 



(135) 



We see that we can compute all such matrix elements for 
an operator A by first computing the object 

J J dxdy(j){z,z,c,d)A(l){z,z,c' ,d'), (B6) 

applying the differential operators -D"^'"~ and D^i^i^^ , 
and then finally setting c,d,c',d' appropriately. Note 
that we have swapped the z and z in <p{z, z,c,d) because 
this generator corresponds to the complex-conjugated 
wavefunction. Since the parameters c and d are always 
real numbers, complex conjugation is implemented sim- 
ply by swapping z and z in the generator ^(z,z). The 
utility of this approach is that we need only compute 
the double integral once and for the simplest possible 
wavefunctions — pure Gaussians. We have outlined the 
approach for single-particle matrix elements, but the 
same procedure also applies to two-particle matrix ele- 
ments such as we have for the Coulomb interactions be- 



tween electrons. 

Each term of the double well Hamiltonian can be com- 
puted in the manner outlined above. Fortunately, the 
integral for the double harmonic potential can be com- 
puted exactly. The integrals over x and y can also be 
computed for the impurity terms once the l/r Coulomb 
potential is rewritten as an exponential using the follow- 
ing identity: 

- = ^^^e-^ (B7) 
r Jo Vs 

This relation also allows the four coordinate integrations 
to be performed in the case of the Coulomb-interaction 

terms. However, for both the Coulomb terms and the im- 
purity terms, the remaining integration over s from the 
above identity must be done numerically for each set of 
quantum numbers since we do not have a closed form ex- 
pression for this integral as a function of the parameters 
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c and d. 

There are several ways that one could treat the kinetic 
terms. In order to compute these terms, we need to first 
act with the kinetic operator on 0(z,z, c, d) as shown 
schematically in Eq. (|B5I) . Since the kinetic operator 
acts only on (j){z , z , c, d) , we can rewrite it as a differen- 
tial operator acting on c and d instead of on z and z. 
Explicitly, it has the form 



IC+ = 



2hiL)4.ijj_ 



huj 



-dddc + hoj^ 



Xo UJ+ 



dc 



d± 



x^h (w+ — 



Xofl 



{uj+-LU-){c-d). (B8) 



the number of matrix elements that need to be computed. 
First of all, each single-particle term in the Hamiltonian 
is real and symmetric: 



(77, n, m|yl|?7', n' , m!) — (?]', n' , m'\A\rj^ n, m) 



(B9) 



We have included the parameter 77 = ±1 in the states to 
designate whether the state is associated with the right 
dot (77 = — !) or the left (77 = -1-1). As far as the impurity 
terms go, this is the only symmetry they obey, at least for 
an impurity located away from a; = 0, 7/ = 0. The kinetic 
and double well potential terms, however, also preserve 
a left-right symmetry which mixes with the single-dot 
parity operator (—1)"+" : 



2. Symmetries {r],n,m\A\r^' ,n' ,m') = {-1)'^+'^' {-7^,n,m\A\-r]' ,n' ,m') . 

(BIO) 

In computing single and two-particle matrix elements. The two-particle Coulomb matrix elements form a sub- 
there are a few symmetries that we may exploit to reduce group of the permutation group: 



(77i,7ii,mi;772,n2,m2|C|77i,7^'i,77i'i;?72,"-2,m2> = (?72,n2,m2;77l,7^l,ml|C|?72,?^2)"7.2;^7l,n^,77^'l) 

= ('7l,'^i,"li;?72,'^2)"7.2|C|77l,7il, 777,1; 772, "-2,^2) = (773, 77-2, 7713; 77^, 7^1, 77li|C|772, ^^2, "12; 771, 77l, 77I1). (Bll) 

The Coulomb matrix elements also exhibit left-right symmetry laced with parity: 

(?7l,r7l,mi; 772,r7,2,TO2|C|7;i,77,i,77li; 772, 772,m2) = (-l)"l+""+"'l+"2 (-771,771,777,1; -772, 77,2, 77l2|C|-7?'i,77i, 777,i; -772, 77,2, 7713) 

, mi 



3. Orthonormal basis 

Once we have the single and two-particle matrix el- 
ements for the various terms in the Hamiltonian, it is 
useful to switch to an orthonormal basis so that we may 
subsequently apply the standard Slater-Condon rules (re- 
viewed in Appendix [C]) to obtain multi-electron matrix 
elements. There must exist a linear transformation which 
takes an operator A in the original basis to A! in the or- 
thonormal basis: 

A! = LAL\ (B13) 

where appears to ensure the hermiticity of A' . Con- 
sider now the overlap operator O: 

{r],n,m\0\ri' ,n' ,m') — {r],n,m\r]' ,n' ,m') . (B14) 

If we apply the transformation to O, then we must obtain 
the identity operator by definition: 

C = LOLt = 1, (B15) 

implying that 

O^L-\L-^)l (B16) 



This factorization of O is called the Cholesky Decom- 
position. The components of L^^ in the original (non- 
orthonormal) basis form a lower triangular matrix, and 
they can be found in a systematic way. The fact that 
is triangular also fits with the fact that it is the transfor- 
mation which takes non-orthonormal states to orthonor- 
mal states. To see this, consider a matrix element of an 
operator A between two non-orthonormal states: 

{N0i\A\N02} = {N0i\L-^LAL^L^)-^\N02) 

^ {N0i\L-^A\L^)-^\N02). (B17) 

From this, it is clear that the orthonormal states are ob- 
tained by applying (U)^^ to non-orthonormal states: 

\0) = {L^)-^\NO). (B18) 

In summary, to switch to the orthonormal basis, 
we first compute the matrix of overlaps in the non- 
orthonormal basis and then perform a Cholesky Decom- 
position on the result to obtain L. For each single- 
particle term A in the Hamiltonian, we compute LAL'^ to 
obtain A' in the orthonormal basis. For the two-particle 
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(Coulomb) terms, we use a slight generalization of the 
transformation to the case of multi-index tensors: 

C'aptS = X! Laa'Li3fj:Ca'li'-f'S'L^^,^L\,^. (B19) 

Here, the index a represents a distinct set of single- 
particle state quantum numbers, a = {77, n,m}, and sim- 
ilarly for /?, etc. 

Appendix C: Slater-Condon rules for multi-particle 
matrix elements 

Once we have all the matrix elements for each term 
in the Hamiltonian in the orthonormal basis, we can ap- 
ply standard Slater-Condon rules to obtain matrix ele- 
ments of the various terms with respect to fully anti- 
symmetrized multi-electron wavefunctions. Since we are 
interested in multi-particle states that have net spin zero, 
we can write a generic multi-electron state as follows: 

|vI.) = |{aW,af',...},{aW,af ,...}), (CI) 

where for example, the symbol a^^^ represents a spin-up 
electron in the Fock-Darwin state with quantum numbers 



QfC^) = {r]k,nk: TTifc}. We define Ivf) to be a fully antisym- 

(k) 

metrized state, so the alj. , Vfc, must be distinct from 

(k) 

each other, and likewise for the a^j^ . Furthermore, the 

particular order of the a^^"^ (and similarly for the a^^^) 
does not matter. Therefore, we may pick a canonical or- 
dering, and for concreteness we choose this ordering to 
be determined by the energy of the single-particle states, 
with lowest energy states to the left (smaller k) (see Eq. 
(lAlSp ). Obviously, states which only differ by the choice 
of rj (which dot they belong to) will have the same energy, 
so we furthermore specify that in these cases, the state 
with 77 = — 1 resides to the left of the state with rj—l. 

We first consider the single-particle terms in the Hamil- 
tonian, which include the kinetic, double well potential, 
and impurity terms. Denoting these terms collectively by 
the operator A, we consider the following matrix element. 



If 1 4*) = l^"'), then according to the Slater-Condon rules, 
we have 



(vI'|^|vI/) = ({a(i\af\...},{4^\af\...}|^|{4^\af\...},{4^Up,...})- ^ {n\A\n) + ^ {k\A\k}. 

kG{q<,",...} Ke{a[^\...} 



If \^>) and 1 5*') differ by only one single-particle state, for example (assuming the states that differ have spin 
either in the set {a^\a^^\ ...} or in then up) 



(vI^VI*) = af\ ...}, {a['\ ...}\A\{a['\ af, ...}, ...}) = (-1)^+^' (af ^ |^|af ), (C4) 

I 



where a^'^ is not in the set {a^^\ a^^ \ ...} and a'j^ 

is not in the set {af\ ...,a\^\ ...}. If |*) and l^-') differ 
by more than one single-particle state, then the matrix 
element vanishes. 



The story is a little more complicated in the case of 
two-particle operators, which for us only includes the 
Coulomb interactions. For the diagonal matrix elements, 
we have 



(^-ICI*) = {{a\'\...},{Jl\...}\C\{a\'\...},{af\...}) = («:A|C|«,A)- J] {k,X\C\X,k)- ^ {^,X\C\X,k) 



K,AG{a<'',...}U{a<'',...} K,Ae{Q<'\...} K,Ae{~<' 



(C5) 
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If \^) and \^') differ by one single-particle state (and assuming this state has spin up), we have 

{^'\m = ({af ' , af\...}, {a«, ...}|C|{a«, «f , ...}, {a«, ...» 

= J2 {K,af^\C\K,af )-{-!)'+'' («,«f'^|C|af (C6) 

Ke{{a\^^\...}n{a'^'> ,...})u{a[^\...} ,...}n{4'\...} 

I 



In the case of two-particle operators like the Coulomb 
interaction, if |^) and \^') differ by two single-particle 

states, then the rcsiilt docs not vanish. Furthermore, 
the result will depend on whether the two single-particle 



states have the same spin or not. We first consider the 
case where they do have the same spin, and suppose for 
concreteness that they both have spin up. Then 



(*'|C|*) = ({af\...,af\...,af\...},{4^\...}|C|{4^...,af\...,af\...},{4^),...}) 



= (_l)^l+^2-|-<-|-4 



a. 



n - (a 



(4) 
t 



(C7) 



A similar relation holds for the case where both pairs of pairs have opposite spin, we instead have 
distinct single-particle states have spin down. When the 



(^'ICI*) 



if'' 



af ^...}|C|{a;^^...,a|-^...},{ai^ 



,-}> 



(-1) 
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